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Abstract 

It is well known that density matrices can be used in quantum mechanics to 
represent the information available to an observer about either a system with a 
random wave function ("statistical mixture") or a system that is entangled with 
another system ("reduced density matrix"). We point out another role, previously 
unnoticed in the literature, that a density matrix can play: it can be the "con- 
ditional density matrix," conditional on the configuration of the environment. A 
precise definition can be given in the context of Bohmian mechanics, whereas or- 
thodox quantum mechanics is too vague to allow a sharp definition, except perhaps 
in special cases. In contrast to statistical and reduced density matrices, forming 
the conditional density matrix involves no averaging. In Bohmian mechanics with 
spin, the conditional density matrix replaces the notion of conditional wave func- 
tion, as the object with the same dynamical significance as the wave function of 
a Bohmian system. 

PACS number: 03.65.Ta (foundations of quantum mechanics) 
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In this paper we shall be concerned with the following claim: Once we deal with 
particles with spin in Bohmian mechanics, we are more or less obliged to regard the 
quantum state of any system (except the universe) as given by a density matrix, which 
then has precisely the same dynamical significance as the wave function. The aim of 
this paper is to elaborate on this statement, as it is far from obvious in what sense a 
density matrix could represent the dynamical state of a Bohmian system. In fact, our 
statement is in sharp contrast with that of Bell j2j: 

So in the de Broglie-Bohm theory a fundamental significance is given to the 
wave function, and it cannot be transferred to the density matrix. 

Although this is correct for spin particles, the situation changes as soon as we consider 
spin or any other internal degree of freedom. To appreciate this point, it is essential to 
distinguish between different roles that density matrices can play in Bohmian mechanics 
(or, for that matter, in other versions of quantum mechanics). In one of these roles, the 
density matrix is of a purely epistemic character, i.e., it expresses ignorance, whereas 
in another role, a role that has as yet not been discussed in the literature and of which 
Bell was obviously not aware, a density matrix is of direct significance to the Bohmian 
particle motion, as the "conditional density matrix." 

We distinguish in this paper five roles of density matrices: the statistical, reduced, 
combined (reduced statistical), conditional, and fundamental density matrix. We ex- 
plain the relations between them and their relevance to the particle motion. We explain 
in particular the new notion of conditional density matrix and its relevance to Bohmian 
mechanics. 

A particular consequence of our discussion is that the same system can, at one and 
the same time, have a conditional density matrix and, say, a different reduced density 
matrix. Thus, when speaking about "the" density matrix of a system, it is necessary 
to specify whether one refers to the reduced or the conditional density matrix. This 
is new: among the traditional types of density matrices, it is always clear (except for 
the ambiguity in some cases as to whether one should consider collapsed or uncoUapsed 
wave functions) which type of density matrix is relevant to a given system, and what 
this density matrix is — so that it is possible to speak of the density matrix of the system. 
The fact that a system can have two different density matrices at the same time is why 
we have to focus on the role that a density matrix plays for the theoretical treatment of 
a system, since that is the only way to understand how more than one density matrix 
can be relevant to the same system. 

2 Bohmian Mechanics 

We begin by briefly recalling Bohmian mechanics. It is a theory of point particles 
moving in physical space M^. For the sake of concreteness, consider a universe of 
nonrelativistic particles whose positions we denote by Qi{t), . . . ,Qj^(t). They move 
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according to Bohm's equation of motion, 



at rrij ip*ip 

where rrij is the mass of particle j, ip : is the wave function, and denotes 

the scalar product in C'^. In the case k = 1 (spin 0), ((T]) simplifies to 

at rrij ip 
ip evolves according to the Schrodinger equation 

j=i 

where the potential V may take values in the k^k Hermitian matrices. The configuration 
Q{t) = . . . , Q^{t)) is random and |?/'(t)p-distributed at every time t, 

FiohiQit) edq) = \tP{q,t)\^dq. (4) 

This is possible because of an equivariance property of (P) and Q: if (jH) holds at 
t = then it also holds at every other time. This follows from the following continuity 
equation, a consequence of 0: 

^ = -divm\) (5) 

where v is the velocity field, i.e., the (time-dependent) vector field on whose j-th 
component is the right hand side of IQ. We remark that the state at time t of a Bohmian 
universe is described by the pair {Q(t),ip(t)). 

What we describe in this paper about conditional density matrices applies not only to 
conventional nonrelativistic Bohmian mechanics as just described, but also to Bohmian 
mechanics on curved manifolds lU]. to Bohm's trajectories for Dirac wave functions 
(see jni p. 272] and |HI), to the photon trajectories of [T7], to the jump processes of [12], 
and, in a sense that we will explain more fully in Section 17.41 also to theories with a 
variable number of particles (TUl [H HH [121 • 



3 Three Density Matrices 

If 7i denotes the Hilbert space of a system S, a density matrix for is a positive, 
(bounded) self-adjoint operator W : Ti ^ 7i with trH^ = 1. If, as in Bohmian mechan- 
ics, 7i is a space of wave functions on a configuration space Q, H = L^{Q,C^), then a 
density matrix can also be viewed as a function W : Qx Q ^ End(C^) (where End(C'^) 
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denotes the space of linear mappings (endomorphisms) C'^ C'^). The translating 
relations between the two views, operator on 7i and function on Q x Q, are 

{W^^Jy{q) = fdq'J2 q') r' (q) and (6a) 

Q 

W:,{q,q') = {q,s\W\q',s') (6b) 
where s and s' index the standard basis of C'^. The function W has the properties 

W{q',q) = W*iq,q') (7a) 
< jdqjdq' ^^:(g) W:,{q, q') r\q') < oo V^^ G H 



Q Q 



(7b) 

dgtrcfc iy(g,g) = 1, (7c) 



where W* denotes the adjoint endomorphism in C*"', whose matrix is the conjugate 
transposed. Conversely, the properties ((Tj) are sufficient for W to define a density matrix 
W . A particular consequence of (f7a|) is that on the diagonal of Q x Q, W{q, q) is a 
Hermitian endomorphism (and thus tr^fc W{q,q) G M), and a particular consequence of 
fl7b|) is that 

trcfe W{q, g) > Vg G Q. (8) 

There are four ways in which density matrices can arise from Bohmian or quantum 
mechanics. Three of them are well known; we briefiy recall them anyway. 

1. First, by statistical mixture. Suppose the wave function ip oi a. system is random 
with probability distribution n{dip) on the unit sphere S^iTi) of the Hilbert space 
Ti. The associated statistical density matrix is 



Vrstat= j /i(rfV')l^)(^| (9a) 
y(H) 

respectively 

w^statMg,?')= / iimr{q)rAq')- w 

This density matrix was first considered in PH] • Note that different distributions fi 
may lead to the same density matrix. (For example, the density matrix on the 
finite-dimensional Hilbert space C'^ arises from the discrete uniform distribution 
over the vectors of any orthonormal basis in C'^, as well as from the continuous 
uniform distribution over the unit sphere ^(C'^).) The significance of VTgtat hes in 
the fact that the distribution of the random outcome Z of an experiment performed 
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on the system depends on /i only trough Wstnu i-e., different /I's leading to the 
same density matrix also lead to the same statistics of outcomes. More precisely, 
when the experiment "measures the observable" A, the probability of obtaining 
an outcome Z in the set i? C M is 

Prob(Z eB)= tr(W^,tatPi(fi)) (10) 

where is the projection-valued measure (PVM) on the real line given by the 
spectral decomposition of the self-adjoint operator A} This follows by averaging, 
according to /i, of the probability that the result is in B given that the state vector 
of the system is ip, which is (in both standard quantum mechanics and Bohmian 
mechanics) {ip\P^{B)\ip) . A particular consequence of (fTUI) is that the outcomes 
of position measurements are distributed according to the density 

p(g) = trcfcVrstat(g,g) (11) 

on configuration space Q. 

From Schrodinger's equation Q for ip, one obtains an evolution law ^18j for pl^stat: 



'^^f^ = [H, W,t.t] (12a) 

respectively 

^^ dW,Uq,q') ^ fy^^^^^^^^^ _ H,,W,,,,{q, q') (12b) 

where Hq means that the Hamiltonian H acts on the variable g, and [ , ] denotes the 
commutator. We remark that VTstat is "pure," i.e., a projection to a 1-dimensional 
subspace, if and only if n is concentrated on that subspace. 

2. The second situation in which a density matrix is relevant involves a system 5*1 that 
is entangled with another system 5*2. In this case, the composite system Si U S2 
possesses a wave function \E''^^*^(gi, ^2) or \E' G TYi 7-^2, but no wave function is 
associated with Si alone. However, the following reduced density matrix can be 
associated with Si. 

W^red = tr2|^)(*| (13a) 

^We remind the reader that in Bohmian mechanics such an experiment need not measure anything 
in the hteral sense of the word 01 E]. We also note that (|10|l holds not only for "measurements 
of observables," but for arbitrary experiments £ with results in the value space V: with every £ is 
associated a positive-operator- valued measure (POVM) IZHH such that the probability of obtaining 
from £ an outcome in the set i? C V is iviW^tatPsiB)) ■ 
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respectively 



W,e4liqu q[) = j dq^Y, ^'^^H^i, ^2) K'.sM^ Q2) (13b) 
Q2 

where tr2 denotes the partial trace over 7^2- This kind of density matrix was first 
considered in jTH]. Note that Wj-ed is an operator on Hi. Like W^stat, W^red possesses 
significance in terms of probability distributions: if one "measures" A on 5*1 alone, 
then the probability of obtaining a result Z in the set S C M is 

Prob(Z G 5) = tr(H^rcdPi(5)) (14) 

where the trace is, of course, taken in Tii. This equation follows from the fact that 
the observable on Tii (g) 7^2 that corresponds to this experiment, as an experiment 
on Si U 5*2, is i (g) i, so that the probability for Z G S is {"^iP^iB) (g) which 
equals (HH). 

If 5*1 and 5*2 are decoupled, i.e., if if = Hi^l + 1^ H2, the reduced density matrix 
evolves in the same way as statistical density matrices do, governed by Hi: 

zh^^ = [Hi,W,,^] (15a) 



respectively 



^^ WedMO ^ Hi,W,,a{q, q') - Hi^'W^q, q') . (15b) 



In case 5*1 and 5*2 are coupled, Wj-ed does not have an autonomous dynamics, i.e., 
its evolution depends on the \l/ from which it arises. We remark that l^red is "pure" 
if and only if Si and S2 are disentangled, \l/'^^''^(gi, ^2) = i'l'iqi) i^2^{q2)- 

3. The third possibility is the combination of the first and the second types of density 
matrices: the reduced density matrix of a statistical mixture. Suppose the wave 
function \l/ of the system Si U 5*2 is random with distribution ^ on yiHi 7^2)- 
Then define the combined density matrix by 

W^comb= J Mc^*)tr2|*)(*| (16a) 

respectively 

vr,o„.b:;(gi,gl)= J /i(rf^) y"rfg2$^^^^^Hgi,g2)^:;.,(g;,g2). (I6b) 
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This kind of density matrix was first considered in [THl p. 424]. PVcomb can be ob- 
tained either by averaging the reduced density matrix associated with the random 
state or by reducing, i.e., taking the partial trace of, the statistical density 
matrix on Tii ® 7^2 associated with fi. Again, the probability that the result Z of 
an experiment on Si "measuring" A lies in the set C R is 

Prob(Z eB) = tr(W^eombPi(5)) . (17) 

This follows either from averaging (fT^ over /i or from applying (fTUj) to A i. 

Like the reduced density matrix, Macomb follows the unitary evolution governed by 
Hi whenever that makes sense, i.e., whenever Si and 5*2 are decoupled. VFcomb is 
pure if and only if fj, is concentrated on the subspace Cipi ® ^2 for some ipi ^'Hi. 



4 A Fourth Density Matrix 

We now turn to the fourth, novel, kind of density matrix: the conditional density matrix. 
It also involves a system Si that is entangled with 5*2, and it is related to the notion 
of conditional wave function which we recall first. For the sake of definiteness, we 
take 5*2 to be the environment of 5*1, i.e., the rest of the universe. 

In Bohmian mechanics for spin particles, more precisely in Bohmian mechanics 
with complex-valued wave functions, the conditional wave function of 5*1 is obtained 
from the wave function \l/(gi,g2) of Si U 5*2 by inserting the actual configuration Q2 of 

V'cond(gi) = ^^(gi,Q2) (18a) 



/AT 

^heie M = j dqi\^>{qi,Q2)\^ (18b) 
Qi 

is a normalizing factor ensuring that / iV^condP = 1- V^cond can be viewed as the wave 
function of 5*1 alone. It does not, in general, evolve according to a Schrodinger equation 
(121), indeed it does not have an autonomous dynamics at all.^ In fact, in appropriate 
situations the evolution of V'cond leads to collapse, in the usual textbook manner, which 
seems quite appropriate for the wave function of a subsystem, ipcond shares the following 
basic properties with the wave function in Bohmian mechanics: 

• The conditional distribution of Qi given Q2 is iV'condP- More precisely, we have 
the following formula for the conditional probability: 

Prob(Qi G dqi\Q2) = |V^cond(gi)Pc?gi , (19) 



^The conditional wave function at time t — need not determine the conditional wave function at 
later times. As an example, consider two situations with the same ^, the same Q2(0) and different 
Qi{0): since V'cond does not depend on Qi, it will be the same in the two situations at t = 0, but since 
the motion of Q2 typically depends on Qi, the two situations will typically have different (52 's at later 
times, and thus typically different V'cond's. 



7 



which resembles the formula (jH) for the probability in terms of the wave function. 
(fTUI) follows from the fact that the pair {Qi,Q2) is \^\'^ distributed. 

• The motion of Qi can be computed from ■i/'cond according to 

^ = ^I.n^i^(Q„.....Q,.J. (20) 

dt mij T/^cond 

which is the same formula as (j21) for the velocity in terms of the wave function. 

An analogous conditional wave function cannot be formed, however, when the par- 
ticles of 5*2 have spin or any other internal degree of freedom entailing that the wave 
function has several complex components. The reason is that ipconA as defined in (I18a|) 
would have too many components, i.e., more spin indices than appropriate for a wave 
function of 5*1 alone. In particular, ipconA would not be an element of Tii. 

We propose to consider instead the conditional density matrix^ which is obtained 
from '^{qi,q2) ^*{.q'iA'2) ^y inserting the actual configuration Q2 of 5*2 for both q2 and 
^2, and contracting over the spin index belonging to S2: 

W^cond:((gi,gO = ^$^^^^^Hgi,Q2)^.V,(g;,Q2) (21) 

with normalizing factor'^ 

M= fdq,J2^^'''iquQ2)^:,JquQ2). (22) 

One easily checks that PVcond satisfies ((H) and thus is a density matrix.^ The expression 
for the corresponding operator VFcond reads 

where Pq^ is the projection- valued measure on Q2 defined by the joint spectral decom- 
position of all position operators of 5*2, and the fraction is a Radon-Nikodym derivative 
of an operator- valued measure on Q2 with respect to a real- valued measure on Q2, and 
thus an operator- valued function on Q2, into which we insert Q2- 

We remark that VFcond is pure if and only if "^{qi, Q2) as an element of L^{QiX''')^ 
C^^ is a tensor product, \l'*^'*^(gi, Q2) = i^i^{qi)i^T- particular, VFcond is pure if \E' is 
complex valued. 

The conditional density matrix has the following properties analogous to those of 
the conditional wave function: 



■^One can show that for almost every configuration Q — (Qi, Q2) (ahiiost every with respect to the 
l^'P distribution), JV will be neither zero nor infinite. 

''The only step that may not be obvious is the finiteness part of (j7b|) . which follows from the fact 
that A/" < 00 so that for any fixed value of S2, '^'^^'^^ {qi, Q2) as a function of si and qi lies in L^(Qi, C^^ ); 
thus the scalar product with any -tp G L^{Qi, C'^'^) is finite. 
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• The conditional distribution of Qi given Q2 can be computed from H^cond by taking 
the trace on the diagonal. More precisely, we have the following formula for the 
conditional probability: 

Prob(Qi G rfgilQa) = tr^fci Wcond{qi, Qi) dqi . (24) 

This follows from the fact that the pair ((5i,<52) is \^\'^ distributed. Note that 
the right hand side is the usual expression (fTT|) for the probability distribution on 
configuration space when a system is described by a density matrix. 

• The motion of Qi can be computed from VFcond according to 

dt mij trcfei PVcond(gi,gi) 
To be able to appreciate , we have to consider a fifth type of density matrix. 



5 A Fifth Density Matrix 

A density matrix is relevant in yet another way: in a modified version of Bohmian 
mechanics in which the particles are guided not by a wave function but by a density 
matrix. Let us call this VT-Bohmian mechanics. Whereas in the conventional version of 
Bohmian mechanics the wave function (of the universe) is something real, as an objective 
component of the state of the universe at a given time, in H^-Bohmian mechanics instead 
of a wave function (of the universe) we may have only a density matrix. This density 
matrix does not arise in any way from an analysis of the theory, but is built into the 
fundamental postulates of W^-Bohmian mechanics. It is a fundamental density matrix, 
Wfund) in contrast to the four other density matrices we have discussed, which were 
derived objects, derived from ip and Q. Like the conditional density matrix, Wfund has 
not been considered previously in the literature. The state at time t of a VT-Bohmian 
universe is given by the pair (Q(t), Wf^^dit)), and it evolves according to 

dQj h Vqtrcfciyfund(g,g'), , 

—f- = — Im ' — (g = q =Q) (26) 

dt rrij tTckW{und{q,q) 

as the equation of motion for Q, and 

^^dm^ = [H, #f„M] (27a) 
ot 

respectively 

^^murAq^ql ^ H^w^^^^^q^ _ H,Wi,ndiq, q') (27b) 
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for Wfund, respectively Wfund(Q'5 Note that equations (jTfj) are the same as (fT^ and 
(lisp . (j26p was first written down by Bell [2^ for the purpose of contrasting it with 
the implications of Bohm's equation of motion (^Q) for a system with a random wave 
function, hence described by VFgtat- 

The configuration Q{t) is random with distribution given by the trace of the diagonal 
of W'fundW, i.e., 

Prob(g(t) G dq) = trcfe ^undlg, q, t) dq. (28) 

This is possible because of the following equivariance theorem: if (I28|) holds at t = 
then it also holds at every other time. To see this, note that (jTTj) implies that 

atrcfe Vrfund(g,g) I. ... I ^ ^ i^^. 

— = - div (trcfc Wiund[q, q) V) (29) 

at 

where v is the velocity field, i.e., the (time-dependent) vector field on Q whose j-th 
component is the right hand side of (j26p . 



6 Discussion 



Bohmian mechanics, as described in Section |2l is a special case of ly-Bohmian mechan- 
ics: if PVfund is pure, i.e., if it arises from a wave function ip via 

Wi^r.6'Aq,q') = riq)rAq'), (30) 

then the equation of motion reduces to Bohm's equation of motion (^, the prob- 
ability law reduces to the \'ip\'^ law (jH), and the evolution (jTTj) entails that PVfund 
remains pure and arises from a wave function that evolves according to the Schrodinger 
equation 0. 

Conversely, the equations (j26p and (j28|) of ly-Bohmian mechanics arise for the be- 
havior of subsystems from Bohmian mechanics for systems of many particles with spin: 
The motion of the particles of subsystem Si is governed according to (j2SI) by a density 
matrix, IVcond, in the same way as in W^-Bohmian mechanics the motion of particles 
is governed according to by a density matrix, IVfund- In addition to the veloci- 
ties, also the probabilities (|24j) are determined by a density matrix in the same way 
as in ly-Bohmian mechanics ()28p. Thus, even were the universe as a whole governed 
by Bohmian mechanics, for most subsystems the state would be described by a density 
matrix, PVcond, with the velocities and probabilities of the subsystem governed by the 
equations of H^-Bohmian mechanics for PVcond- In this sense, l^-Bohmian mechanics is 
the theory relevant to most systems in a Bohmian universe. (More precisely, this holds 
for all those systems for which Vl^cond is not pure.) 

A big difference, however, between the dynamics of a subsystem and VT-Bohmian 
mechanics lies in the fact that, unlike the fundamental density matrix, see (j27j), the 
conditional density matrix need not evolve unitarily. Nevertheless, there are special 
situations in which VFcond does evolve unitarily, at least as a good approximation. This 
happens trivially when Si and 5*2 are disentangled, \E'(gi,g2) = i'liqi) ® 'ip2{q2), and 
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decoupled (so that they stay disentangled). It also happens when (and for as long as) 
5*1 and 5*2 are decoupled and 

^{qu q2) = Mqi) ® ^2(^2) + ^^(gi, ^2), (3i) 

i.e., 

^^^^-{q,,q,) = ^r(gi)^f (gs) + (^^)^^^Hgi, ^2), (32) 

where iIj2 and have disjoint g2-supports and Q2 G support ^/'2- Such a situation often 
occurs after a measurement, and indeed allows us to regard ipi as the (effective) wave 
function of Si, obeying Schrodinger's equation Q. For spin 0, (jH^ characterizes the 
situation in which we can expect the conditional wave function to evolve unitarily; thus, 
the conditional density matrix evolves unitarily in all situations in which the conditional 
wave function would for spin 0. We obtain another case of unitarily evolving IVcond by 
replacing by 

^'^'%qi,q2)=ri'''{qi)Mq2) + {^^r''{quq2), (33) 

with a complex- valued ^2, and assuming in addition that the Hamiltonian H2 for S2 
involves no interaction between spin and configurational degrees of freedom. 

For example, consider an EPR-Bohm-Bell pair of spin 1/2 particles, each headed 
towards its Stern-Gerlach magnet, with qi and q2 the positions of the particles. Suppose 
both magnets are oriented so as to measure and that the geometry is such that particle 
1 completely passes its SG magnet before particle 2 reaches its SG magnet. Initially the 
spin state is the singlet state, depending on neither qi nor q2, and we may assume as well 
that the configuration space wave packet is initially of product form ipi{qi)ip2{q2) ■ Then 
the initial wave function is of the form (j33|) with \I'-'- = and (regarding the possible 
values of Sj as ±1) 

ri^'^iqi) = -^{Ss,,iSs,,-i - (5,„_i(5,,,i)Vi(gi) , (34) 

corresponding to 

^(IT)li)-|i)IT))®V'i (35) 

in the standard representation. 

Until particle 1 reaches its magnet the Schrodinger evolution preserves this form and 
W^ond = ® where ipi = ipiif) obeys Schrodinger's equation for particle 1. 

Moreover, until particles 2 reaches its magnet (i.e., in the absence of a magnetic field 
acting on particle 2), H2 involves no coupling between spin and translational degrees 
of freedom, so that the form (p?!^ is preserved and IVcond evolves unitarily according to 
P7jl . even after particle 1 has reached its magnet. After particle 1 has passed trough its 
magnet (but before particle 2 reaches its magnet) Wcond = |(W^up + Wdown), where W^p, 
respectively Wdown, corresponds to the pure state | T ) ® '^/'up, respectively | i ) ® V'down, 
the states to which || ) ® t/^i, respectively | i ) ® V^i, would evolve under the Schrodinger 
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evolution for particle 1. After particle 2 reaches its magnet, W^cond no longer evolves 
unitarily (or even autonomously). Rather it collapses either to W^p or Wdown according 
to whether the initial configuration is such that Q2 ends up going down or up. 

Throughout the course of the entire experiment Qi evolves according to (j^^ . (Note 
also that after particle 2 has crossed its magnet, is again approximately satisfied, 
with the wave packet that does not contain Q2.) 

We now turn to the relations between the various density matrices, and discuss first 
the relation between VFcond and W^ed- ^red is the average conditional density matrix, 
with the average taken with respect to quantum equilibrium, i.e., over the ensemble in 
which Q = {Qi,Q2) is l^&P distributed: 

W,edl\{qi,q[) = j c?gic^g2 1^(^1,^2)1' W^eond:((gi,gl)(g2). (36) 
S1XQ2 

This relation makes clear that a system can have a conditional and a reduced density 
matrix at the same time, the two being different from each other: the conditional 
density matrix of a system depends on the configuration Q2 of its environment; when 
this dependence is averaged out by taking the quantum equilibrium expected value one 
obtains the reduced density matrix of the system. (Note that for spin (jHBj) is the 
quantum equilibrium average of |'0cond)('0cond|-) 

Similarly, the combined (reduced statistical) density matrix is an average of the 
conditional density matrix, with the average taken over the ensemble in which \& is /i 
distributed and, given ^f, Q is distributed: 

Macomb:; (gi,g;)= j ^^{d^>) j rfQi^Qs |^(Qi,Q2)|' w^cond:;(gi,gl)(Q2,^). (37) 

y{H) Si X Qa 

Of course, Wstax can also be viewed as an average (of over the ensemble with 

/i-distributed but this does not involve the conditional density matrix. 

The fact that VTcond determines the Bohmian velocities according to should be 
contrasted with the failure of such a connection for VTstat, W^red, and Wcomb'- If the wave 
function ^ of a system is random, the Bohmian velocities have to be computed from the 
actual reahzation of tp, and thus could assume different values, corresponding to different 
T/^'s, even when Q is held fixed. Inserting, for example, W^stat in a formula like ()25|) or 
fl2^ would yield, in contrast, an average velocity at Q, averaged over the ensemble of 
different ip's (with the additional Q-dependent weight proportional to \il){Q)\'^). This is 
what Bell referred to in the phrase we quoted in the beginning, and what he elucidated in 
|2] . Similarly, since W^ed is the average of the conditional density matrix, over a certain 
ensemble, it leads to an average velocity (in fact to the best guess at the velocity that 
one could make without knowing Q2)- In contrast, PVcond depends on the actual value 
of Q2 and yields the true Bohmian velocity, as defined by and the wave function of 
the universe. 

The statistical analysis of Bohmian mechanics in remains valid when conditional 
wave functions are replaced by conditional density matrices. 
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7 Remarks 



7.1 Conditional Density Matrix in Orthodox Quantum Me- 
chanics 

In orthodox quantum mechanics, the definition (|2ip of the conditional density matrix 
cannot be written down, for lack of a configuration Q2 that could be inserted into \E'. 
However, orthodox quantum mechanics arguably maintains that macroscopic objects 
can be viewed and treated classically, which presumably means that there should exist 
something like a "macroscopic configuration." In case that ^ is such that the conditional 
density matrix does not change much with the microscopic details of Q2 (i.e., that it 
is quite accurately determined by merely the macroscopic information about (52), a 
conditional density matrix also makes sense in orthodox quantum mechanics. In this 
case the conditional density matrix of orthodox quantum mechanics would equal, within 
its accuracy, the one of Bohmian mechanics. Another way of obtaining this density 
matrix is to collapse the wave function (to the region of configuration space having 
q2 compatible with the actual macroscopic configuration of 6*2), and then to take the 
reduced density matrix. 

7.2 Second Quantization 

In ^21, we describe a construction that might be called the "second quantization of a 
Markov process." Parallel to the "second quantization" algorithm of forming a Fock 
space out of a given 1-particle Hilbert space and the free Hamiltonian on Fock space 
out of a given 1-particle Hamiltonian, this construction builds a dynamics on the con- 
figuration space of a variable number of particles out of a given 1-particle dynamics. A 
key step in this construction is a general procedure for forming the law of motion for 
particles, given an arbitrary 1-particle law. Interestingly, the conditional density matrix 
is indispensable for this procedure (except when wave functions are complex-valued). 

A Bohm-type law of motion for one particle associates a velocity vector field on 
with every (smooth) 1-particle wave function. We now regard this association abstractly 
as a given mapping, from which we want to systematically construct the A^-particle law 
that provides the velocities of all particles from an A^-particle wave function and the 
positions of all particles. By inserting the positions of all but one particle into the wave 
function, we get a conditional object for one particle — for spin a conditional wave 
function, otherwise a conditional density matrix. Only if the one-particle law associates 
with this conditional object a velocity vector field on M^, can we insert the position of 
the remaining particle into the vector field and get the particle's velocity. For spin > 
we thus need more than what we mentioned at the beginning of this paragraph: we need 
that the one-particle law provide a velocity field for every density matrix, as ly-Bohmian 
mechanics does, and not merely for every wave function. 
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7.3 Empirical Consequences of VK-Bohmian Mechanics 

One may wonder whether one can decide empirically between Bohmian mechanics and 
Ty-Bohmian mechanics, or, in other words, whether one can determine empirically in a 
universe governed by VT-Bohmian mechanics if the fundamental density matrix is pure 
()30|). The question is delicate. We think that the answer is no, for the following reason: 
compare a VT-Bohmian universe with a Bohmian universe with a random wave func- 
tion such that the associated statistical density matrix equals the fundamental density 
matrix of the IV-Bohmian universe. Since an empirical decision, if it can be made at 
time to, would have to be based solely on the configuration Qt^ at that time, and since 
the distribution of Qt^ is the same in both situations, it seems that there cannot be a 
detectable difference: A given Qt^ could as well have arisen from an appropriate wave 
function from the random wave function ensemble as from the corresponding fundamen- 
tal density matrix. 

What makes the question delicate, however, is, in part, the following: we might not 
take seriously a theory involving a wave function of the universe or a density matrix 
of the universe that is "unreasonable" or "conspiratorial." Therefore, the question is 
connected to questions such as what would count as a "reasonable" VFfund, and whether 
a statistical mixture mimicking a given "reasonable" VFfund might have to contain some 
"unreasonable" wave functions. 

7.4 Conditioning on Spatial Regions 

It is often desirable to define the subsystems Si, i = 1,2, as encompassing all those 
particles which are presently located in the regions Ri C M^, with i?i U i?2 = and 
Ri n R2 = 0. To condition on the configuration Q2 of 5*2 then means to condition on 
the configuration in the region R2. We describe below what appears to be the most 
convenient way to carry out such a conditioning on a spatial region. One might suspect 
that conditioning on a spatial region is a very complicated story. But, in fact, it could 
not be simpler. 

Since the number of particles in the region Ri can vary over time, it is helpful to 
consider right from the start a configuration space of a variable number of particles. We 
consider the space 

00 

r(M3) ;= IJ R^^/Sn (38) 

ra=0 

where Sn denotes the group of permutations of n objects, which acts on M^" by permuting 
the particle labels. A configuration from r(]R^) represents any number of identical 
(unlabeled) particles. For a discussion of this space, see [T^ . 
We can extend the definition of F to arbitrary sets R, 

00 

r{R) := U R^/Sr,. (39) 

n=0 
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When i? C M^, T{R) can be viewed as a subset of r(M^), containing those configurations 
for which all particles are located in R. Now observe that, when Ri and R2 are disjoint 
sets, then 

r(i?i U R2) = T{Ri) X r(i?2) . (40) 

This property is helpful, as it tells us that the definition of the subsystems Si in terms 
of spatial regions Ri leads to a Cartesian product decomposition Q = Qi x Q2 of 
configuration space, and thus allows us to use, without change, all of our considerations 
on conditional density matrices, which assumed such a decomposition. 

8 Conclusions 

We have introduced the notion of conditional density matrix in Bohmian mechanics, 
and contrasted it, on the one hand, with the notion of conditional wave function, and 
on the other hand, with various other notions of density matrices. In contrast to the 
statistical, reduced, or combined (reduced statistical) density matrix, the conditional 
density matrix possesses direct significance for the particle velocities. 

The fact that with the same system can be associated several density matrices brings 
into sharp focus that the meaning of a density matrix is not a priori; instead, various 
meanings are conceivable. Ultimately, the meaning of a density matrix arises from its 
relevance to the primitive objects, such as particle world lines, that the theory is about. 
In Bohmian mechanics, the various types of density matrices that we have considered 
are all relevant to the particles, but in very different ways. 
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